arXiv: 1506.07195v2 [physics.plasm-ph] 4 Aug 2015 


Under consideration for publication in J. Plasma Phys. 


1 


The Wisconsin Plasma Astrophysics 
Laboratory 

C.B. Forest^t, K. Flanagan^, M. Brookhart^, M. Clark^, C.M. 
Cooper^, V. Desangles^, J. Egedal^, D. Endrizzi^, I. V. Khalzov^, H. 
Li^, M. Miesch®, J. Milhone^, M. Nornberg^, J. Olson^, E. Peterson^, 
F. Roesler^, A. Schekochihin®’^, O. Schmitz^, R. Siller^, A. 
Spitkovsky®, A. Stemo^, J. Wallace^, D. Weisberg^ and E. Zweibel^ 

^Department of Physics, University of Wisconsin-Madison, Madison, WI 53706 USA 

^Laboratoire de Physique de I’Ecole Normale Superieure de Lyon, CNRS and Universite de 
Lyon, 46 allee d’ltalie, 69364 Lyon Cedex 7, France 

^National Research Centre “Kurchatov Institute”, Moscow, 123182, Russia 
^Los Alamos National Laboratory, Los Alamos, NM 87545, USA 

®High Altitude Observatory, National Center for Atmospheric Research, Boulder, CO 

80307-3000, USA 

®Rudolph Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX16NP, UK 

^Merton College, Oxford OS14JD, UK 

^Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08540, USA 
(Received xx; revised xx; accepted xx) 

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility 
designed to study a range of astrophysically relevant plasma processes as well as novel 
geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed 
from strong samarium cobalt permanent magnets now confines a 10 m^, fully ionized, 
magnetic-field-free plasma in a spherical geometry. Plasma parameters of Tg « 5 to 20 
eV and rig « 10^^ to 5 x 10^^ cm“^ provide an ideal testbed for a range of astrophysical 
experiments including self-exciting dynamos, collisionless magnetic reconnection, jet 
stability, stellar winds, and more. This article describes the capabilities of WiPAL, along 
with several experiments, in both operating and planning stages, that illustrate the range 
of possibilities for future users. 


1. Introduction 

During the past five years, a medium-scale multi-investigator experimental 
plasma facility—the Wisconsin Plasma Astrophysics Laboratory (WiPAL)—has been 
constructed and is now in operation at the University of Wisconsin-Madison. At the 
heart of the facility is a high bay that houses a 3 m diameter spherical multi-cusp 
confinement device (see figure [^. Its flexibility and diagnostic access allow WiPAL to 
currently support two major experiments: the Madison Plasma Dynamo Experiment 
(MPDX) and the Terrestrial Reconnection Experiment (TREX), along with several 
new experiments being added in the near future. Both the flexible design of the facility 
and the shared hardware between experiments allow for a rapid turn-around between 
different configurations. 
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Figure 1: The WiPAL facility (left) consists of a main laboratory space where the vessel is 
housed and several auxiliary spaces for high voltage power management, water pumping, 
computer control, and housing the TREX insert. Photo in upper right shows the 3 m 
diameter vessel covered with arrays of probes and clear viewports. 


The purpose of this paper is to describe 
the physical parameter regimes that can 
be achieved in WiPAL, to illustrate its 
flexibility, and to show how the plasma 
confinement scheme lends itself to many 
experiments of interest to the astrophysics 
and fundamental plasma physics communi¬ 
ties. Already, this facility has transitioned 
from a single-purpose experiment focused 
on plasma dynamos (MPDX) to a multi¬ 
investigator facility accommodating a new 
reconnection experiment (TREX). In the 
future, WiPAL will transition into a 
collaborative user facility open to other 
laboratory astrophysics research groups. 
This paper illustrates the potential of the 
facility by describing (in addition to MPDX 
and TREX) several additional experiments 
that are now being pursued: acoustic and 
Alfven wave propagation in connection 
with helioseismology; pulsar and stellar 
wind launching from a rotating dipolar 
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Figure 2: A static calculation of the vac¬ 
uum magnetic field from the permanent 
magnets at the edge and the background 
earth field. This calculation has been 
normalized by hall probe measurements 
along a radial chord. 
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high-power helicity injection. In addition to these experimental efforts, numerical and 
theoretical computations, often carried out with identical dimensionless parameters, are 
used to create predictive models and inform elements of design. 

This article begins with an overview of WiPAL along with experimentally achieved 
plasma parameters. Then a description is given of the diagnostic suite and of the magnetic 
cusp confinement used in WiPAL. This is followed by brief summaries of experiments 
that are under way or in development. 


2. Description of the Facility 

The WiPAL facility (shown in figure consists of the plasma confinement vessel and 
the associated infrastructure (high bay, electrical power, water cooling, and associated 
laboratories). The vacuum vessel consists of two 1.5 m radius cast-aluminum hemispheres 
mounted on a track that can be separated during a vacuum opening. This provides the 
opportunity for quickly inserting different devices inside the WiPAL vessel. For example, 
a transition from full MPDX operation to the TREX configuration requires only a few 
days. 

Nearly 3000 samarium cobalt (SmCo) permanent magnets (each with |B| >3 kG 
at their surface) are held in axisymmetric rings on the inside of the vessel. These rings 
alternate in polarity to form a strong high-order multipole field that decreases in strength 
to the background earth field within 20 cm of the vessel wall (shown in figure . This 
edge-localized cusp field provides sufficient plasma confinement to achieve the parameters 
listed in tablej^while leaving the ions in the core unmagnetized. The rings were carefully 
produced to assure axisymmetry, thus reducing many aspects of WiPAL experiments to 
two dimensions. Since all the magnetic field variation is in the radial and polar directions, 
gradient and curvature drifts are azimuthal. This leads to an azimuthal symmetry of the 
generated plasmas, thereby reducing convective losses associated with magnetic ripple. 

To date, argon and helium plasmas 
have been created and heated by 
an array of hot emissive lanthanum 
hexaboride (LaBg) cathodes, each 
independently powered by separate 
30 kW power supplies (currently 12, 
with 18 planned). See figurefor an 
example discharge in a typical helium 
plasma. Five 20 kW magnetrons at 
2.45 GHz are in the process of being 
added to independently heat electrons 
through electron cyclotron resonance 
(EGR) heating in the magnetized 
edge, bringing the maximum steady- 
state input power to ~ 650 kW. These 
sources are quasi-stationary, meaning 
that plasmas can be sustained for 
tens of seconds limited only by vessel 
heating. 

A unique capability of the facility 
is the ability to control the plasma 
rotation at the boundary (Gooper 
et al. 20141. Early experiments 
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Table 1: Table of parameters achieved to date 
in separate discharges (with 10 of 18 possible 
cathodes and no EGR heating) as well as 
projected parameters. The ionization fraction, 
/%, is calculated using neutral pressure at the 
vessel wall and line-integrated electron density. 
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Figure 3: Time traces of plasma parameters and input power from a typical WiPAL 
discharge in helium (left) and argon (right). Radiated power is measured with a 
bolometer, is measured with a pressure gauge at the edge of the vessel (assuming 
room temperature neutrals), Tg is measured with a Langmuir probe inserted into the 
unmagnetized bulk, and rie is measured with the mm-wave interferometer system. 


highly conducting, flow-dominated plasmas and understanding the interface between 
magnetized and unmagnetized plasma conditions. Plasmas have been stirred at speeds up 
to 10 km s“^ in the azimuthal direction by applying a J x B torque in the magnetized edge 
region. Each of the cathodes are mounted on a motorized stage so their insertion depth 
can be quickly and accurately set between shots. Cathodes pushed into the unmagnetized 
bulk of the plasma ionize and heat the plasma, while cathodes pulled back into the 
magnetized region draw current across magnetic field lines, injecting torque into the 
plasma. This torque is viscously coupled to the unmagnetized core. This boundary-driven 
flow scheme was designed for dynamo studies on MPDX, but will be used to study other 
features of fast-moving plasmas in the WiPAL facility. 

Additionally, several capacitor banks with ignitron switching have been used to drive 
short pulse reconnection experiments in TREX ((3.31 and to inject current for plasma 
jet experiments ((3.6). Eor example, initial reconnection experiments utilized an up to 
10 kV, 100 fjF capacitor bank to drive approximately 20 kA for 250 /xs, generating a 10 
G reconnection field from a two-turn internal coil. WiPAL is also equipped with a large 
Helmholtz coil pair which provides a uniform held throughout the plasma volume of up 
to 275 G. 


2.1. Major Diagnostics 

The WiPAL vessel is covered with approximately 200 ports for diagnostic access of 
plasma discharges via robotic probe arrays and advanced optical measurements. Owing 
to the high temperature and steady-state nature of WiPAL discharges, in situ probes 
have been carefully designed to accommodate high heat hux. Additionally, the small loss 
area provided by the cusp conhnement is increased signihcantly by inserting probes into 
the plasma. Motivated by these issues, a suite of non-invasive optical diagnostics has 
been developed and works in concert with inserted probes. 

2.1.1. Two-dimensional Probe Systems 

Most plasma astrophysics experiments require high-resolution, two-dimensional (2D) 
measurements obtained using internal probe arrays. To meet this need, WiPAL has 
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several robotically controlled probe drive systems in place. Two types of scanning systems 
are used, as shown in figure The first probe system is a set of sweeping probes (one 
currently installed, with six planned) that are inserted through a single ball-joint vacuum 
seal port, combining radial and angular motion to sweep out a 2D plane. Two stepper- 
motor-controlled transverse stages allow the regions shown to be scanned shot-to-shot 
with arbitrary resolution. The second probe system is an array of single-axis probes 
that are inserted every 10° in latitude. These probes have a single stepper-motor drive, 
allowing radial shot-to-shot scans. 

The main probe type to be used in 
WiPAL is a combination velocity-magnetic 
probe. Each probe tip consists of two Mach 
probes (four planar Langmuir faces collecting 
ion saturation current) to measure two 
orthogonal components of v. Development 
is under way to add three Hall sensors to 
measure the entire three-dimensional (3D) 

B field in each probe tip. These probes are 
designed to withstand long, warm plasma 
pulses, utilizing a thermally conductive 
copper shaft paired with an electrically 
and thermally insulating quartz shield. The 
lifetime of probes is ultimately limited by 
accumulated heat load over many plasma 
pulses, but thermocouple temperature moni¬ 
toring facilitates optimization of plasma duty 
cycle versus probe temperature. 

Increased probe coverage also leads to 
diminished plasma performance because 
plasma losses to the exposed probe area can 
be much larger than the edge losses through 
the permanent magnet cusp field. Motorized 
probes allow for flexibility in controlling the plasma loss area. For example, staggered 
radial scans of various probes can keep the loss area constant while still mapping out the 
full extent of probe coverage. 

The 2D structure of the magnetic field on the surface of the plasma will be measured 
using an array of 64 external magnetic probes placed in the region between the plasma 
and the vessel wall. These probes will consist of three-axis Hall probes for measuring 
low-frequency magnetic fields as well as two orthogonal Mirnov coils for measuring 
higher-frequency fluctuations. The final implementation of this full array will resolve 
£=8 harmonics in the polar direction and m=8 harmonics in the toroidal direction. 

In addition to this collection of internal probes, WiPAL is outfitted with an array of 
optical diagnostics. A millimeter-wave interferometer system measures the line-integrated 
electron density; a Fabry-Perot spectrometer provides measurements of ion temperature; 
and optical emission spectroscopy (OES) can be used to determine neutral and ion density 
profiles, electron temperature, and electron density. All of these systems provide precise, 
non-invasive measurements of key plasma parameters without increasing the plasma loss 



Figure 4: An example of a predicted 
dynamo-capable axisymmetric flow, with 
poloidal flow (left) and toroidal flow 
(right) is well diagnosed by the full 
motorized probe suite. WiPAL probe 
coverage will include both sweep probes 
(6 total, left) and linear probes (18 total, 
right). 


area. 
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2.1.2. Absolute Density Measurements 

WiPAL is equipped with a heterodyne millimeter-wave (320 GHz) interferometer 
system designed to measure chord-integrated absolute electron density via wave phase 
shifts. This system is very similar to interferometers used on the Helically Symmet¬ 
ric Experiment (HSX) and the Madison Symmetric Torus (MST) ( |Deng et al. |200^ 
20061.Time resolution of density measurements is set by the intermediate frequency (IF) 
of the variable-frequency source. Typically, the IF is set to / ^ 1 MHz, but is adjustable 
from 0.1-100 MHz. 

At WiPAL densities of rie = 10^^ to 5 x 10^^ cm“^, phase shifts of several fringes (one 
fringe corresponds to a 27r phase shift between the reference and plasma beams) must be 
measured. Additionally, owing to time resolution requirements, this measurement must 
be made at frequencies up to several megahertz. This constraint is met by using a high¬ 
speed held-programmable gate array (FPGA) programmed to compute the phase at each 
period of the reference beam. The ultimate density resolution is set by the speed of the 
FPGA clock relative to the IF frequency and yields Sn ^ 10^° cm“^ for an IF of 1 MHz. 

This absolute electron density measurement is used to calibrate Langmuir probe data 
for point measurements as well as to constrain neutral emission models for spectroscopic 
measurements of the electron temperature ( ^2.1.4 1. Because of the good time and den¬ 
sity resolution, this diagnostic will be used in TREX to measure density fluctuations 
associated with magnetic reconnection. 

2.1.3. Ion Temperature and Flow 

Doppler-shifted line emission is measured using a Fabry-Perot spectrometer. Plasma 
light is passively collected from either the A = 488nm argon ion line or the A = 468.6nm 
helium ion line complex depending on the working gas. By imaging these ion lines at 
a high resolution, the Fabry-Perot captures the line-integrated ion velocity distribution 
function. Ion temperature is inferred from the thermal broadening of the lines, while 
velocity can be measured via the Doppler shift of the peak away from its stationary 
value. 

The Fabry-Perot system has several distinct advantages over a grating-type spectrom¬ 
eter. For a comparable sized grating, the Fabry-Perot has an increase of nearly 100 
times the resolving power. Additionally, because the output of the Fabry-Perot is a 
symmetric ring structure where cx A, integration around the rings can greatly increase 
the signal-to-noise ratio. This procedure, called ring summing, results in a sensitivity gain 
of 10-30 compared to linear cuts of the pattern (Goakley et al. [1996 ). Additionally, it 
allows for relatively short integration times (roughly 1 s in this system) using a standard 
digital camera. WiPAL also has a high-performance intensified charge-coupled device 
(CCD) camera capable of taking even shorter exposures (approximately 0.1s) for multiple 
measurements during discharges. The typical resolution (i.e., error bars) of Fabry-Perot 
ion temperature measurements is STi « 0.05 eV and the velocity resolution is of the 
order <5?; « 10 m s“^. 


2.1.4. Electron Temperature and Distribution 

Several survey spectrometers are used to measure emission throughout the visible 
spectrum with lower spectral resolution than the Fabry-Perot. The spectrometers capture 
segments of the near-ultraviolet to the near-infrared range, either monitoring the entire 
spectrum at low resolution (300 nm — 888 nm at 0.387 nm) or small segments at a higher 
resolution (381 nm — 511 nm at 0.045 nm). Spectra are continuously captured through 
WiPAL discharges with exposure times of 2-80 ms and rates of 10-50 Hz which provide 
a time history relevant for long-pulse time scales. 
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The OES system, calibrated by the millimeter-wave absolute density measurement, can 
estimate the electron temperature and ionization fraction. Line ratios are compared to 
collisional radiative models calibrated to data from other experiments and used to infer 
Te while including additional corrections due to non-Maxwellian electron distribution 


functions in argon plasmas 1 

Boffard et aZ.||2010 

Wang et aL|2013 

I. In helium, line ratios 

can also be used to determine Tg and rig with similar modeling ( 

Schmitz et al. 2008 

). Good 


agreement has been found between the models and the data in the regimes calibrated 
for tokamak scrape-off layers (high Tg, rie). Work is under way to benchmark these 
models for more relevant parameters. In addition to using the OES system for its own 
measurements, WiPAL offers an opportunity to extend these calibrations to regimes 
relevant for laboratory plasma astrophysics. 


3. Major Experiments in WiPAL 

3.1. Large Multi-Cusp Plasma Confinement 

WiPAL is the largest axisymmetric magnetic ring cusp ever constructed for confining 
plasmas. The permanent magnets effectively limit the loss area to ^ 1% of the total 
vacuum vessel surface area, which in turn leads to high, steady-state confinement. 
The combination of WiPAL’s cusp confinement and diagnostic capabilities presents an 
opportunity for studying and characterizing the confinement of magnetic cusps. The cusp 
field is created by 36 rings of alternating-polarity magnets, as shown in figure This 
field is localized to the edge, dropping off to the background earth field within 20 cm 
of the vessel wall. Plasma losses in this configuration are limited to a small cusp width 
on the face of each magnet. To electrically isolate the plasma from the vessel wall, each 
magnet is covered by a thin insulating ceramic tile. 

The long established empirical loss width for a magnetic cusp is the hybrid gyroradius 
w ~ y/PiPe (Hershkowitz et al. |I975| ), yet newer studies have found a scaling with neutral 


pressure which suggests more detailed physics than was previously employed (Hubble 


et al. 20141. Since WiPAL is able to create and sustain plasmas at very low pressures, this 


neutral pressure scaling can be investigated over several orders of magnitude. Accurate 
measurements of the ion velocity distribution at the cusp can be made with the Eabry- 
Perot system. Additionally, owing to the insulating boundary condition, in situ probes 
at the ceramic tile can directly measure the loss width. Understanding the loss width 
scalings in a magnetic cusp is important not only for predicting plasma confinement in 


WiPAL but also to a number of other applications such as plasma processing (Malik 


et al. 

1994 

1979 



and neutral beam injectors (Stirling et al. 


3.2. Self-Exeiting Dynamos 

Astrophysical plasmas are often characterized by high magnetic Reynolds numbers 
{Run = VL/t], where rj is the magnetic diffusivity, U is a characteristic flow speed, and 
L is the system size) wherein turbulent, flow-dominated plasmas continuously transform 
kinetic flow energy into magnetic energy. Understanding this energy transformation and 
predicting what form the magnetic field might take, be it small-scale turbulent magnetic 
fields or large-scale magnetic flux, is the dynamo problem. The theory of dynamos has 
shown that the scales at which the magnetic energy grows are largely determined by 
the relative values of the magnetic Reynolds number and the fluid Reynolds number 
{Re = VLjv, where v is the dynamic viscosity). The ratio Pm = Rm/Re, called the 
magnetic Prandtl number, is thought to be a critical parameter that governs the nature 
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Figure 5: a) Hydrodynamic simulations of edge-driven, two-vortex axisymmetric flow 
in He for increasing Re. Axis of rotation and symmetry is vertical; left hemispheres 
show poloidal flow stream lines and poloidal flow magnitude, right hemispheres show 
toroidal flow contours. Poloidal flow scale is amplified 4x relative to toroidal flow scale, 
b) Kinematic dynamo growth rate of the flows shown in (a). Increasing Re results in 
a larger amount of kinetic energy in poloidal flows, which in turn lowers the critical 
Rm required for positive magnetic eigenmode growth. Flow is only axisymmetric for 
Re < 300; hydrodynamic instabilities produce non-axisymmetric modes at higher Re, 
requiring full 3-D solutions to the induction equation. 


of many astrophysical phenomena since it sets the relative collisional dissipation scales 
of the fields and flows. Astrophysical plasmas span a wide range of Pm; diffuse plasmas 
have Pm ^ 1 whereas denser plasmas have Pm ^ 1. 

Dynamos can be classified as small-scale or large-scale. Small-scale dynamos tend 
to generate magnetic energy but little net magnetic flux, whereas large-scale dynamos 
generate both net flux and energy. While the process by which turbulence generates 


magnetic energy at small scales seems theoretically well understood ( Schekochihin et al. 
2004 Iskakov et aL|[2007 ), understanding how a large-scale magnetic field self-organizes 


from small-scale magnetic fluctuations in astrophysical systems remains a grand challenge 
for plasma astrophysics. 

Studying plasma dynamos in the laboratory requires a previously unexplored regime 
of laboratory plasmas. Unmagnetized, fast-flowing, and highly conducting plasmas are 
required so that magnetic fields can be stretched and amplified by the plasma inertia. Cast 
in dimensionless terms, dynamo action requires plasmas with high Rm and high Alfven 
Mach number {Ma = VIVa ^ 1). The particular geometries being pursued in the MPDX 
configuration of WiPAL build upon similar geometries used by liquid-metal experiments. 
Mechanically stirred liquid-sodium experiments have observed spontaneous magnetic 
field generation ( Monchaux et al.]|2007 2009 Gallet et al.||2012 | and have added to our 
understanding of astrophysical and geophysical dynamos (Lathrop & Forest 2011). A 
plasma experiment has the potential to extend these studies to parameters more relevant 
to astrophysics. Beyond the obvious fact that most naturally occurring dynamos are 
plasmas, the use of plasma rather than liquid metals corresponds to magnetic Reynolds 
numbers increased by a factor of 10 or larger. Additionally, viscosity can be varied 
independently of the conductivity, with Pm ranging from 0.1 to 10, reflecting the wide 
range found in different astrophysical systems. 
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Figure 6: Numerical solutions for turbulence driven small scale dynamo and evidence for 
large scale field generation. To achieve Rm = 2000 and Re = 400 corresponds to a helium 
plasma with V^max=10 km s“^, Te=30 eV, Ti=1.2 eV, rig = 10^® m“®, and =1.2. 
The left plot shows where simulations find turbulent dynamos in Re and Rm space. 



Extensive modeling and theoretical work have been carried out to find flow schemes 


for MPDX that excite a large-scale dynamo (Spence et al. 2009 
2012a|6 


Khalzov et al. 2013 


Katz et al. 2012). By biasing the cathodes in the magnetized edge, torque is 


injected at the boundary of the plasma and is viscously coupled to the core. Dynamo¬ 
relevant flows are found by solving the Navier-Stokes equation for a particular boundary 
condition imposed by the cathode drive. The resulting flow is then used to solve the 
kinematic dynamo equation. One such flow that results in dynamo action is shown in 
figure The imposed boundary condition drives two counter-rotating vortices in each 
hemisphere with a small flow direction reversal near the equator. For fairly modest Rm 
and Re, positive dynamo growth is expected. Work on MPDX is currently directed at 
optimizing the flow drive and attempting to create this two-vortex flow. 

Fast small-scale dynamos can be excited by chaotic flows with the high values of Rm 
and Pm possible in WiPAL. Chaotic flows with positive Lyapunov exponents could be 
achieved with large Re turbulence or by using time-dependent and highly viscous laminar 


flows (Khalzov et al. 2013). Three-dimensional numerical simulations using turbulent 


flows where Re ~ 500 and Rm is very high have confirmed that turbulent, small- 
scale magnetic fields naturally develop in high-Pm plasmas. Moreover, the magnetic 
amplification of these flows is roughly independent of Rm as Rm continues to increase, 
a hallmark of fast dynamos. Figure shows several example of flows where fast dynamo 
growth is possible. Simulations predict fast dynamo growth at plasma parameters that 
can be achieved with increased input power from ECR heating. Because Rm is very 
large in most astrophysical systems, virtually any astrophysical dynamo must be fast. 
Understanding the conditions for fast dynamos and exploring whether they can generate 
magnetic fields on large scales is an open problem in plasma astrophysics. 


3.3. Magnetic Reconnection 

One of the fundamental processes in nearly all magnetized plasmas is magnetic re¬ 
connection: a change in the magnetic topology of the plasma that converts stored 


magnetic energy into particle kinetic energy (Priest & Forbes 2000). Given its importance. 


magnetic reconnection has been studied extensively through theory, numerical computa- 
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Figure 7: Toroidal field coils and a cylindrical insert holding additional coils for driving 
reconnection (left) will be inserted into the WiPAL vessel. The large set of external 
Helmholtz coils seen in the figure (right) is a recent upgrade to the WiPAL facility with 
utility for many future experiments. 




Figure 8: Reconnection phase space diagrams for weak guide-field reconnection (left) and 
strong guide-field reconnection (right) identifying various regimes of reconnection. The 
star characterizes the first plasma obtained with TREX, demonstrating its ability to 
access the new collisionless reconnection regime. 
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tion, spacecraft observations, and laboratory experiments. In recent years, new frontiers 
have emerged with an emphasis on particle energization, reconnection with multiple 
X-lines in large systems, the role of kinetic effects in collisionless reconnection, and 
extending the evolution of reconnection from two dimensions to three. For experiments 
to remain relevant and to advance this maturing field, new devices are needed that access 
these regimes. To address this need, we are now operating the Terrestrial Reconnection 
Experiment (TREX), which is the largest dedicated reconnection experiment to date. In 
addition to WiPAL’s existing vacuum vessel, plasma confinement, heating, and diagnostic 
suite, TREX consists of a cylindrical insert holding coils used for driving reconnection 
(shown in figure]^. The parameter regimes expected for TREX are shown in table 
The full operation of TREX within WiPAL renders the facility unique in its ability to 
address the expanding frontiers of reconnection research. 

An important component to the implementation of TREX into the WiPAL user 
facility is the cylindrical insert housing the reconnection drive coils and the addition 
of a toroidal field coil. The insert and TF coil allow TREX to access multiple magnetic 
configurations (e.g. antiparallel reconnection, strong guide-field reconnection, and 3D 
reconnection), as well as a wider range of reconnection regimes. TREX also utilizes a 
suite of probe arrays specially developed to characterize reconnection dynamics. These 
include stationary magnetic and electric probes along with sweep probe arrays that allow 
for a 2D region to be mapped out over repeated shots ( ^2.1.1 ). The necessary dynamics 
can be captured on time scales up to 20 MHz with the diagnostic suite. For a desired 
value of the reconnecting field, the configuration must be driven at the corresponding 
loop voltage VJoop = 2TrR{0.1vABr). For example, a loop voltage of 5 kV for a 150 G 
reconnecting field is required to access fast reconnection in a hydrogen plasma. This 
consideration has necessitated the development of reliable, high-power pulsing circuits 
utilizing class D ignitron switches to drive reconnection. These drive circuits can be 
pulsed as fast as every 10 seconds. 

So far, the term “collisionless reconnection” has referred to systems where the electron 
and ion distributions can remain near Maxwellian, but the collisionality is sufficiently low 
that their collective fluid behaviors decouple at the ion scale and Hall currents become 
important. However, in a truly collisionless plasma, pressure anisotropy develops which 
strongly impacts the properties of the reconnection process in ways not accounted for 


in traditional Hall reconnection. In fact, spacecraft observations ( 

Hwang et al. 

2013) 

and kinetic simulations at the full ion-to-electron mass ratio ( 

Le et al. 

2015 

show 


that large-scale current layers are driven by electron pressure anisotropy that builds 
in the reconnection region due to kinetic electron trapping effects. To maintain pressure 
anisotropy, the time between electron collisions must be long compared with the full 


transit time of a fluid element through the reconnection layer (Egedal et al. 2013 


Le 


As a valuable tool for displaying the various regimes of reconnection and their tran¬ 
sitions, Daughton & Roytershteyn (2012) developed the reconnection phase diagram 
spanned by the Lundquist number S and the normalized system size A with respect to 
the ion sound Larmor radius, ps = \/mi{Te + Ti)/eB, or the ion skin depth, di = cjiUpi, 
for strong and weak guide-field reconnection, respectively (Ji & Daughton 2011). A 
convenient way of representing the constraint for anisotropic pressure on a system is 
the condition S > lQ{mi/mf.){L/di). The anisotropic pressure region of this phase space 
is shown in figure]^ As indicated by the star, TREX has already demonstrated its ability 
to access this regime of collisionless reconnection. TREX is able to experimentally study 
the role that electron pressure anisotropy has on particle heating. In addition, the narrow 
current layers driven by the anisotropy may in 3D be unstable to reconnection at oblique 
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ne(10i® m-3) I Te (eV) | Br (T) | Bg (T) | L (m) 


Terrestrial Reconnection Experiment 
(TREX; UW Madison) 

0.1- 10 

8-40 

0.04 

0-0.3 

0.8- 1.8 

Magnetic Reconnection Experiment 
(MRX; Princeton) 

2-100 

5-10 

0.03 

0-0.1 

0.3-0.8 

Versatile Toroidal Facility 
(VTF; MIT) 

0.1 - 1 

8-30 

0.01 

0.1 

0.3 


Table 2: Key parameters for various reconnection experiments in hydrogen plasmas 
where B^. and Bg are the reconnecting and guide fields, respectively. 



Figure 9: A simulation of solar-like flow inside the WiPAL vessel with streamlines of 
poloidal flow (left). The internal rotation profile of the Sun inferred from helioseismic 
inversions ( Thompson|2014 1. Pink/red denote fast rotation and blue/green denote slower 
rotation as indicated by the color bar. Inversions near the rotation axis are unreliable 
and are thus omitted from the plot. The image of the angular velocity profile is reflected 
about the rotation axis in order to illustrate the full spherical geometry. 


angles. This effect of 3D reconnection may be important to the self-consistent evolution 
and generation of reconnection with multiple X-lines. 


3.4. Acoustic Waves, Helioseismology, and Angular Momentum Transport 

The WiPAL facility can study both basic solar plasma phenomena as well as long¬ 
standing mysteries explored by the heliophysics community, such as angular momentum 
transport in the Sun and other late-type stars. It is well established from observations 
and modeling that late-type stars with convective envelopes spin down as they age due 


to torques exerted by magnetized stellar winds (Mestel 1968 Barnes 2007 Matt et al. 


20151. It is less well understood how this angular momentum loss is transmitted through 


their convection zones and into their radiative cores. Helioseismic probing of our own Sun 
suggests that this coupling is efficient in the sense that the radiative core currently has 
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a rotation rate that is comparable to the convective envelope (Thompson et al. 2003). 
However, these same helioseismic inversions reveal that the radiative core is rotating 
nearly uniformly while the convection zone rotates differentially, with a ~ 30% decrease 
in angular velocity from equator to pole ( Thompson et aZ.|2003 |. The transition between 
these two distinct rotation regimes is called the solar tachocline, which overlaps with the 
base of the convection zone and the convective overshoot region. 

Angular momentum transport within the solar convection zone is attributed to tur¬ 
bulent magnetized convection and occurs on convective time scales of weeks to months. 
By contrast, the transport of angular momentum across the solar tachocline, coupling 
the core and envelope, is thought to occur on much longer time scales of millions to 


billions of years (Spiegel & Zahn 19921. Our current understanding of both short-term 


and long-term angular momentum transport in the Sun relies on incomplete and indirect 
observational data (e.g. helioseismic inversions) and numerical magnetohydrodynamic 
(MHD) models. WiPAL offers an opportunity to study some aspects of the relevant 
dynamics in the laboratory for the first time. 

Using the drive system as MPDX, a solar-like flow can be programmed at the boundary 
of the WiPAL vessel. A calculation of the structure of this flow is shown on the left 
of figure This flow is aimed at mimicking the helioseismic inversion, shown on the 
right of figure with strong toroidal rotation near the equator. Careful scrutiny of how 
such a boundary-forced differential rotation spreads into the interior through meridional 
flows, viscous diffusion, thermal gradients, and magnetic tension will provide insight into 
core-envelope coupling and the resulting spin-down of late-type stars. Magnetic tension 
in particular is thought to be responsible for maintaining the uniform rotation of the 


radiative interior by suppressing shear, as described by Ferraro’s theorem (MacGregor 


& Charbonneau 1999). Furthermore, inserting a dipole electromagnet into the core of 
the WiPAL vessel will create a mock tachocline, testing the structure and stability of 


the boundary layer proposed by magnetic tachocline confinement models (e.g., Gough & 
McIntyre||1998D . 


By exciting the acoustic normal mode spectrum of a spherical plasma with external 
antennas, the stationary plasma response can be characterized. It is predicted by theory 
and observed in the Sun that these normal modes split when rotation is introduced 


(Christensen-Dalsgaard 2002). Helioseismology relies on these mode splittings at different 
locations on the Sun to discern the flow profiles of the interior. A similar analysis is 
available at WiPAL with the addition of the external magnetic probe array described in 
§2.1.1[ When differential rotation is imposed at the boundary, the resulting spectra can be 
spatially decomposed and compared to the stationary mode structure. A mathematical 
inversion can then be used to propose likely candidates for internal global flow profiles, 
which can be benchmarked by Mach probe measurements under similar conditions. This 
would form the basis of a minimally invasive, global flow profile diagnostic at WiPAL 
akin to helioseismic inversions. 

Using this helioseismic diagnostic and solar flow drive, WiPAL will study the short¬ 
term angular momentum transport in the solar convection zone. Experimental set-ups, 
such as that shown in figure can be used to investigate the physics of turbulent trans¬ 
port. In the Sun, turbulent angular momentum transport is thought to be responsible not 
only for sustaining the solar differential rotation, but also for regulating the amplitude 


and structure of the meridional flow by means of gyroscopic pumping (Miesch & Hindman 


2011). More generally, turbulent angular momentum transport, mean flows, and thermal 


gradients in the solar convection zone are all thought to be intimately linked through 
nonlinear feedbacks that can be explored with WiPAL. 
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Figure 10: (a) WiPAL facility with centrifugal wind experiment installed and (b) proposed 
configuration showing the electromagnetic stirring scheme and Y-point. 


3.5. Centrifugally Driven Stellar Winds 

Plasma wind launched from star surfaces carries stellar magnetic field into the local 
heliosphere and out to the interstellar medium. This system is complex and is governed 
by a range of processes, including magnetic reconnection, turbulence and particle heating, 
all topics of present-day heliospheric research. One aspect that has received little experi¬ 
mental attention is the magnetic topology of the advected magnetic field in the interface 
region between the magnetically dominated corona and the flow-dominated wind. This 
region is critical both in creating the Parker spiral and in determining the mass loss rates 
of stars. The Centrifugal Wind Experiment (CWE) at WiPAL will explore centrifugal 
breakout of wind from spinning dipolar magnetospheres. 

The winds of rapidly rotating giant stars and pulsars are particularly relevant. For 
rapidly rotating giant stars, centrifugal breakout has been predicted to proceed in episodic 
bursts of plasmoids ( ud Doula et aL|[2006 |. Recent observations by the Microvariability 
and Oscillations of STars (MOST) telescope have called into question the validity of this 
breakout model, suggesting instead a smooth outflow ( Townsend et qZ.|[ 2013 ). Modeling 
centrifugal breakout in a laboratory setting can advance our understanding of this 
phenomenon. Pulsars have been extensively modeled using analytical MHD and 3D 
particle-in-cell (PIC) simulations of both aligned and obliquely rotating magnetospheres 


(Spitkovsky 2006 Philippov & Spitkovsky 2014). Both cases predict the creation of a 


magnetic Y-point where the closed magnetic topology breaks into an open configuration. 
This is expected to occur at the Alfven radius, the location where the kinetic and mag¬ 
netic energy densities are equal. WiPAL’s large volume, well-developed power supplies, 
and comprehensive diagnostics make it suitable for experimental analysis of these active 
areas of research. 

A laboratory analog of a rotating, stellar magnetosphere is possible using the electro¬ 
static stirring techniques developed for MPDX with the addition of a dipole magnetic 
field at the center of the vessel. Through ECR heating, the dipole magnetosphere will 
be filled with confined plasma. Biased electrodes at the equatorial surface of the dipole 
will establish a cross-field potential gradient, leading to E x B stirring. As the plasma 
spins to critical velocity, a centrifugal wind will be launched. In this way, we can use a 
large cross-field potential to spin the magnetosphere and model the above astrophysical 
situations. 

In preliminary experiments, an electrically insulated, spherical SmCo magnet (r ~ 10 
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cm, surface field ^ 4 kG) will serve as the dipole source, shown in figure]^ ECR heating 
will produce plasma on the spherically symmetric 875 G surface. Using helium, expected 
temperatures and densities will be similar to prior MPDX experiments (Tg « 5 — 20eU 
and He « 10^^ cm“^). Equatorial cathodes, spaced « 10 cm apart, will be biased up to 
1 kV to establish a cross-field potential, producing axially aligned E x B rotation. 

WiPAL’s existing diagnostics can be used for investigation of the Y-Point and current 
sheet regions. The standard suite of Langmuir and Mach probes will be used for basic 
plasma measurements at the edge and will supplement the non-invasive measurements. A 
2D, motor-controlled, 20 MHz resolution magnetic probe array will serve as the primary 
diagnostic. The magnetic probe will look for the characteristic Y-point opening of the 
magnetic field lines at the Alfven radius and will explore the equatorial current sheet at 
larger radii. 

Planned upgrades will explore other regimes. For example, obliquely rotating pulsars 


generate larger spin-down power than aligned rotators (Spitkovsky 2006) and provide 


a more realistic setting for studying current sheet reconnection. By supplementing the 
dipole magnet with phased coils, an obliquely rotating magnetosphere can be created. 
This will allow for experimental confirmation of the above prediction and characterization 
of the change in the undulating current sheets. 

3.6. Stability of Astrophysical Jets in Background Plasma 
Astrophysical jets are collimated magnetized plasma outflows from accreting bodies 
such as active galactic nuclei, binary systems, and young stellar objects. It is thonght 
that jets could play a key role in angular momentum transport in accretion systems. 
Jets are launched and collimated in regions much smaller than the current resolution 
limits of observations, so the details of this process remain elusive. Gertain models of 
the launching process predict that astrophysical jets arise due to magnetic fields in the 
accretion disk ( Lovelace|1976 Blandford fc Payne|1982 Spruit|[2M0 ). In these models, a 
dipole-like magnetic field is sheared by Keplerian rotation in the disk. This creates a large 
electric field between the center and edge of the disk, driving current along the dipole 
field. Then J x B forces accelerate and collimate plasma, creating the jet structure. These 
jets should be susceptible to current-driven instabilities, but recently published theories 
suggest that the presence of external plasma pressure promotes collimation and enhances 
the stability of the jet structure out to large distances (|Lynden-Bell||2003 Li et aZ.|2001 
20061 ). _ 

creates 


The Permanent Magnet Pulsed Jet Source (PM-PJS), shown in figure 11 


an astrophysical jet-like plasma utilizing a design similar to that used by Hsu and 


Bellan (Hsu & Bellan 20021. A large voltage potential is created between two coplanar 


annuli (approximately 30 cm in diameter). This simulates the potential created via the 
aforementioned accretion disk shearing action. The background magnetic field in the 
astrophysical system is simulated using two antisymmetric rings of permanent magnets. 
The use of permanent magnets instead of electromagnets allows the source to maintain 
WiPAL’s cusp-confining boundary condition. Neutral gas is injected via internal channels 
in both annuli with ports placed at eight azimuthally symmetric pairs of magnetic field 
line foot points. This gas is ionized by the large voltage creating eight filamentary loop 
structures. These loops merge and collimate, driving a magnetic jet into the center of 
WiPAL. 


This experiment complements Galtech experiments (Hsu & Bellan 2002 2003 Moser 


& Bellan 2012) by studying the evolution of a magnetic jet evolving into a vacuum. 


Additionally, WiPAL allows the creation of a magnetic jet that propagates into a high- 
/3 background plasma. This will help clarify the relationship between jet stability and 
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Figure 11: Experiment to study Jet expansion into background, finite pressure plasmas. 
Simulation of Jet is taken from Li et al. (2006) and superimposed into WiPAL plasma 
schematic, showing coaxial injection geometry. Right: photo of the assembled pulsed jet 
source preparing for first tests. 


external plasma conditions (Lynden-Bell 2003 Li et al. 2001 2006). Future work will 


focus on the shock and precursors that are formed between the magnetized jet and 


the unmagnetized background plasma as predicted in Li et al. { 2006) and on how the 


development of the kink instability could depend on the strength of the background 
pressure. Experimentally it will be possible to vary the target plasma pressure by more 
than two orders of magnitude. The properties of the expanding magnetic plume (speed, 
density, magnetic field) as well as the shock precursors will be measured using scanning 
probes similar to those used in the reconnection and centrifugal wind experiments. 


3.7. Helicity Injection and Decaying Turbulence 
The generation of large-scale magnetic fields is a fundamental feature of many as- 
trophysical systems. This dynamo action (as described in ^3.2[ ) is often attributed to a 
turbulent upscale transfer of small-scale magnetic helicity, H = J A-BdV. The upscale 
transfer process is ideally described as a local, self-similar inverse cascade where helicity 
is conserved and transferred to larger scales (Frisch et al. | [T^ [ji| [T^ [Blackman fc 
Ji 2006). In real systems however, the transfer process is much more complex and can 


involve non-local transfer directly from the small to large scales as described in the tur¬ 
bulent alpha effect. The upscale spectral transfer of magnetic helicity has been explored 


extensively via direct numerical simulations (e.g. Alexakis et al. 2006). Understanding 
this upscale transfer of magnetic helicity is key to explaining the creation of the large- 
scale magnetic fields observed throughout the Universe. 

In addition to upscale spectral transfer of helicity, transport in space is necessary 
for large-scale dynamo action. Large-scale dynamos can be sustained in systems with 
boundary conditions allowing the outflow of helicity. Roughly speaking, ejection of one 
sign of magnetic helicity allows helicity of the opposite sign to grow without violat¬ 
ing helicity conservation. This idea was originally suggested on theoretical grounds by 


Blackman & Field (2001) and Vishniac & Cho (2001). Later, it was shown that magnetic 


eruptions associated with solar activity transmit a net helicity flux, suggesting that the 
solar dynamo meets these boundary conditions (Rust|2002 Kusano|2002 Liu et al.|2014 


Pevtsov et aZ.|[20l4). In addition to the solar dynamo, this boundary condition effect is 


thought to occur in galactic dynamos driven by supernovae (Rafikov & Kulsrud 2000). 
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There is evidence from numerical simulations that the generation of large-scale fields is 


promoted by boundary conditions that permit escape of magnetic helicity (Brandenburg 


& Sandin 

2004 

Kapyla et al. 2010 

Hubbard & Brandenburg 

2012 


ideas have been tested in the laboratory. 

In the plasma physics community, helicity injection is primarily explored as a startup 


mechanism for magnetically dominated tokamak plasmas (Raman et at. 2003 2010 


Battaglia et aL|2009 |. Some laboratory astrophysical experiments that inject net helicity 


into plasma systems have also focused on magnetic reconnection and other instabilities 


that occur during Taylor relaxation of low-/3 plasmas (Cothran et al. 2009 Gray et al. 
2010[ Jarboe 2005). All of these experiments rely on helicity injection at large scales 


(comparable to the scale of the plasma volume) relaxing via magnetic instability to 
minimum-energy, helicity-conserving Taylor states in low-/3 plasmas ( Taylor|[l986 |. As a 
complement to these studies, we will inject helicity at small scales and use the large size 
of the WiPAL vessel to directly observe the upscale spectral transfer and transport of 
helicity in a high-/?, dynamo-relevant plasma conditions. 

Generalizing the PM-PJS experiment discussed 
above ()3.6), helicity can be injected using the edge- 
localized cusp field inside WiPAL. A large potential 
can be applied between alternating rings, which will 
drive a current along cusp field lines and lead to 
an injection of helicity into the system without a 
guide field present (figure 12) . By using all of the 
rings of permanent magnets inside the vessel, this 
helicity will be injected at a small (^=18, m = 0) 
scale relative to the size of the vessel. The resulting 
rings of flux will then be blown off the wall into 
the bulk unmagnetized plasma with kinetic and 
magnetic energies near equipartition. In order to 
drive an upscale transfer instead of simple resistive 
diffusion, we require the resistive decay time to be 
much longer than the Alfven time of this system. 

Gast in dimensionless terms, this means that the 
Lundquist number, S = TtIta, must be large. Using 
the high-power capacitor banks created for TREX, 
magnetic fields of R ~ 100 G can be induced from 
current driven along the cusp field lines. For WiPAL 
discharge parameters, this corresponds to S' ~ 1000. 

Under these conditions, we expect helicity to undergo 



Figure 12: The axisymmetric ring 
cusp lends itself to a novel he¬ 
licity injection experiment where 
currents, flowing between the 
magnetic pole faces can be 
configured to inject small scale 
magnetic helicity. 


a turbulent transfer before resistive diffusion can dissipate the injected magnetic energy. 
Observing the transfer and transport of helicity in this set-up will complement dynamo 
studies conducted by MPDX and provide the astrophysics community with a laboratory 
environment to probe this fundamental process. 


4. Conclusion 

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) provides the plasma astro¬ 
physics and fundamental plasma physics communities with a unique opportunity to study 
plasma phenomena in a laboratory setting. Hot, dense, unmagnetized, and fully ionized 
plasmas are routinely created and confined in quiescent states for seconds as astrophysics 
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experiments are performed. This user facility has been designed and operated with the 
goal of maximizing both plasma performance and flexibility of use. 

To date, WiPAL has confined steady-state plasmas which are hotter and denser than 
any other large-scale non-magnetically dominated plasma device. These parameters have 
proven to be sufficient for both studying dynamo relevant regimes (MPDX) and providing 
high-Lundquist-number target plasmas for reconnection studies (TREX). In the longer 
term, the heating power will more than double as the cathode system is completed and 
the ECR system is installed. 

Diagnosing the high-performance unmagnetized plasmas in WiPAL has been an area 
of intense focus. Arrays of robotically controlled magnetic and electrostatic probes are 
capable of mapping out large areas of WiPAL plasmas. Advanced optical diagnostics 
have been developed with the goal of reducing the additional loss area added by this 
suite of in situ probes. Millimeter-wave technology is used to power a compact heterodyne 
interferometer system for measuring electron density with high resolution. A Fabry-Perot 
interferometer is used to extract the ion velocity distribution. Finally, a complement 
of low-cost spectrometers coupled with collaborative modeling makes reliable estimates 
of the electron temperature and ionization fraction. All of these diagnostics provide a 
valuable set of data for every WiPAL discharge regardless of experiment. 

Perhaps most significantly, the WiPAL facility has already demonstrated the abil¬ 
ity to quickly change experimental configurations. This paper has outlined two major 
experiments already running on WiPAL (MPDX and TREX) as well as several new 
experiments in various stages of planning and implementation. We envision WiPAL 
transitioning in the near future to a user facility model in which investigators from 
outside the collaboration groups listed above could apply for and receive support to 
carry out experiments with technical support from the WiPAL staff. 

The construction of the facility was supported by a National Science Foundation (NSF) 
Major Research Instrumentation grant. The MPDX and TREX research is now supported 
by NSF and the U.S. Department of Energy (DoE) and the NSF Center for Magnetic Self 
Organization in Laboratory and Astrophysical Plasmas (CMSO). The helioseismology 
studies are part of a collaboration with NCAR (M. Miesch) and are supported by a NASA 
graduate fellowship (E. Peterson). The stellar wind experiment is a collaboration with 
Princeton University (A. Spitkovsky) and is supported by an NSF graduate fellowship 
(D. Endrizzi). 
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